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ABSTRACT 

Optical  properties  of  photonic  crystal  (PC)  structures  have  attracted  much  attention  in  the  field  of  research  area 
because  of  its  ability  to  direct  and  control  light  propagation  in  compact  device  with  feature  sizes  comparable  to  the 
wavelength  of  light.  Therefore,  from  the  early  days  of  quantum  mechanics,  physicists  have  tried  to  understand  the  optical 
properties  of  photonic  crystal  in  Hilbert  space  using  Fourier  expansion. 

In  this  paper  a theoretical  study  of  two-dimensional  (2D)  planar  optical  photonic  crystals  is  carried  out;  this 
optical  photonic  crystal  consisting  of  a dielectric  layer  and  dispersion  less.  This  software  was  written  in  MATLAB  to 
simulate  and  analysis  the  photonic  bands  for  2DPC  of  TEM  polarizations.  Fourier  coefficients  for  the  expansion  of 
dielectric  constant  are  calculated  with  refractive  index  ,7]  = i andra,  = 3.45 . 

KEYWORDS;  2D  Photonic  Crystal,  Fourier  expansion,  Hilbert  space,  Photonic  Band  Gap 

INTRODUCTION 

The  structures  of  photonic  crystals  (PCs)  are  constants  with  dielectric  periodically  modulated  whose  distribution 
follows  a periodicity  of  the  order  of  a fraction  of  the  optical  wavelength  [1].  If  electric  permittivity  depends  on  two  spatial 
variables  only,  ,rand  y,  the  PCto  be  two-dimensional  (2D)  [2,3].  PC  is  a material  that  is  periodically  modulated  the  index 
of  refractive  on  a length  scale  comparable  to  the  desired  operation  wavelength.  It  is  named  a crystal  because  the  crystal  is 
formed  by  a periodic  arrangement  of  basic  building  blocks.  In  addition,  the  term  photonic  is  added  since  photonic  crystals 
are  designed  to  affect  the  properties  propagation  of  photons. 

When  wave  propagation  is  enters  into  a material,  it  made  most  important  affected,because  some  feature  of  interest  of  this 
wave  is  modulated.  The  coherent  waves  are  scattering  at  the  interfaces  between  different  featured  areas.  In  the  case  of 
photons  this  feature  is  index  of  refractive.  Other  systems  propagating  such  as  waves,  sound  or  electrons  within  a 
semiconductor,  these  features  are  the  Young’s  module  and  the  electric  potential  respectively  [4,5]. 

The  optical  photonic  crystals  have  been  the  subject  of  intense  investigation  because  of  its  ability  to  control  the 
properties  of  photons  [6,7].  Photonic  crystals  ability  has  led  to  mold  the  flow  of  light  resulted  in  a variety  wonderful 
fascinating  optical  phenomena,  such  as  omnidirectional  reflection,  low  loss  bends,  high-Q  cavities,  negative  refraction,  and 
the  design  of  thermal  emission  [8]. 

Figure  lshow  the  basic  theoretical  background  of  PC. Lord  Rayleighwas  first  time  studied  the  propagation  of 
electromagnetic  wave  in  periodic  media  in  1887.  This  is  consistent  with  the  one-dimensional  (ID)  PC,  and  he  knew  the 
fact  that  they  have  a narrow  band  gap  prohibiting  light  propagation  through  the  planes,  whichdependent  on  the  angle. 
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because  of  the  different  periodicities  experienced  by  light  propagating  at  cases  of  non-normal  incidences,  resulting  a 
reflected  color  that  differ  sharply  with  the  angle.  Although,  during  the  next  century,  multilayer  films  received  more 
intensive  study,  in  1987,  Yablonovitch  and  John  joined  the  tools  of  classical  electromagnetism  and  solid-state  physics,  that 
the  concepts  of  omnidirectional  photonic  band  gaps  in  2D  and  3D  was  introduced  [9,10].  The  generalization  inspired  the 
name  Photonic  Crystal  was  led  to  many  developments  in  fabrication,  theory,  and  application,  from  integrated  optics  to 
negative  refraction  to  optical  fibers  that  guide  light  in  air  [9,1 1]. 

1-D  2-0  3-D 

mm 

pwgfcln  ^pwte<»cln^ 

Figurel:One,  two  and  three  direction  of  photonic  crystals  periodic.  The  crystal  material  is  periodicity  structure 

[9,10]. 

2d  Photonic  Crystal  Tem-Polarized  Light 

The  electric  field  vector  wave  equation's  for  2DPCsshown  below  [12]: 

Vx  Vx  E(r)  = k2£r(r)E(r)  ^ 

£(r)  = £(x,y)  = £0£r(x,y ) 

E{r)  = E(x,y,z) 

where  £ ( t y)isPC  relative  permittivity  function's.  The  vector  of  the  electric  field  is  a function  of  three  spatial 
coordinates  system.  For  a TM-polarized  field  i .e.'.E(r)  = E (r)z  > the  divergence  of  the  field  is  always  zero.  In  the  normal 
propagation  case,  i.e.  (angle  of  90°  with;  axis),  after  that,  the  electric  field,  /T  ( jq , is  shown  as  a function  of  x and  v,  i.e. 
[12]: 

V 2£.  {x,  y)  + k2 £r  (x,  y)E,  (x,  y)  = 0 ^ ) 


Using  variablesseparation  method,  electric  field  can  be  written  as  below  [6]: 

^rZ(x)  + ^r[X(x)  + ]z(x)  = 0 

!W(y)  + ^[T(: v)-^V(y)  = o 

dy~  c~ 

Equation  (2)  reduces  to  two  separate  equations. One  2D  structure  has  been  divided  into  two  separate  ID  structures, 
one  in  the  vertical  and  the  other  in  the  horizontal  direction. 


Such  as  above  procedure  to  reduce  the  wave  equation  into  two  separate  differential  equations  on  the  lengthof  each 
coordinate,  and  we  are  trying  to  do  this  for  magnetic  TE-polarized  field.  Wave  equation  for  two  dimensions  TE- 
polarization  can  be  considered  below  [6]: 


Vx 


£r(r) 


= k^H(r)  (4) 


£(r)  = £(x,y)  = £0£r(x,  y) 
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e (x,y)  isPC  relative  permittivity  function's.  Vector  of  the  magnetic  field  is  a function  of  the  three  spatial 
coordinates  system.  For  TE  polarized  light  only  the  normal  component  of  the  field  would  be  nonzero.  Below,  the  derive 
ofmagnetic  field  component  equation[6]: 


V rjxiV  xH)  -T]X2H  = k2H(r) 


'?  = 


1 


...(5) 


er(r) 


3 3 


3 3 


+ +TJ^Hz  +*tth' +f?^rHt=-kiH 

dx  ox  dy  dy  dx  dy  dz 


...(6) 


Using  the  Fourier  transformation  for  solution,  so,  the  solution  write  as,  h _{x,y,z)  = A(x,y)e~‘n  where  y coefficient 
corresponds  to  the  wave  vector  angle's  with  respect  to  z axis  [6]: 

3 3 , 3 3 . 32  32  . „ , . 2,  ...  (7) 

— rl  — A + — r}  — A + r}  — A + T]  — A = Mri]-k;) 
dx  dx  dy  dy  dx  dy 


Assuming*  and  y are  field  components,  the  product  of  two  distinct  functions  of  x and  y,  become,  Aix,  y)  = %{x)i//( y) 


and  [6]: 


i a2  , , i a'  „ nli)  n 1 i i ,,  ■■•(8) 

rrZ (■*)  + T^rrYh ’)+  -Tr>7 l -^7  [ z-i/s(y)]  = (r 

In  3v  J V(y)  3y 


ZM  3x2 


v(y)  dy2  \jl  ^ J z(x)  dx 

If  the  wave  of  incident  propagates  (angle=  90 °z  axis  ( y = ()  ))  in  the  plane  of  PC,  and  the  permittivity  to  be  a sum  of 
two  terms,  each  being  single  functions  of  .rand  y,  the  magnetic  field  will  take  the  form  in  equation  (9).  Thus,the  set  of  TE- 
polarization  wave  equation  in  a two  dimension  photonic  crystal  for  in-plane  propagation  found  as  [6]: 

1 3 


dx1  ^ X (a)  + ¥ ( y ) dx  |_  dx 


3 / ^ 

x-X(x) 


Kny) 1 — y 

dy2  X(x)  + Y(y)  dy 


dy 


+ [i;iW-Ak(j)  = 0...(9) 
+[k02Y(y)  + &]y(y)  = 0 


HILBERT  SPACE  FOURIER  EXPANSION  METHOD 

Now,  in  the  case  of  TM  modes  in  2D  PC,  the  study  has  the  two  dimensionsHelmholtz  equation  [2]: 


3 V 3 VI 


y-+yr\-^-y)r 


...(10) 


Wherespectral  parameter  ,j  and  index  of  refractive  nix)  is  a continuous  positive  function  [2]: 
n(x  + mlal+m2a2)  = n(x),  reft",  mt,m2eZ  (H) 

Seek  (//as  a solutions  bounded  nontrivial  of  the  Helmholtz  (equation  (10))  which  is  the  Laplacian  distribution x2y/ 
is  also  borderedin/U . Therefore,  the  solutions  take  Bloch  representation  [2]: 

iy(x)  = eitx</>(x)  (12) 


where,  <p  satisfies  the  periodicity  condition: 
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</>(x + m,a,  + m2a2 ) = fax),  xe  R2,  m, ,m2eZ,  and  the  vector  k = ^ ^ for  certain  integers  f] , ^ and  ^ being  the 
reciprocal  vectors  basis.  Substituting  (12)  into  (10)  we  get:  -V2fax)-2ik.'Vfax)+\\k(fax)  = i)n(x)2fax)  ■ 

After  construct  an  orthonormal  basis  of  the  complex  Hilbert  space  consisting  of  those  functions  in  l2(A)  that 
satisfy  the  above  condition.Writing:#^^  +;i£i  :(tltt2)e z2),  the  orthonormal  basis  is{m(A)"1,2^(  :te  Z2}>  where: 
<p,(x)  = = (/|,/,)eZ!,  t»(A)  = det  A • 

Each  basis  function  satisfies  the  periodicity  condition:  <pt(x  + mlal+m2a2)  = <p,(x),  (m„n»2)e  Z2,te  Z2 . as  well  as  the 

conjugation  symmetry:  = ^(x),  /e  Z2  • 


The  computed  of  the  expansion  coefficients  ^ in  fax)— \ ftp  (x)  shown  as  follows: 

reZ2 

$t  = m(Ay'  J dxfax)(p,  ( x ) = /n(A)-1  J dxfax)(p_, (x)  ’ 

A A 

Observe  that:  _ v>,  (x)  = \\tfy  + t2b2  f<p,  (x)  = A7t2tT  (ATA)~'  up,  (x)  - where  t is  the  column  vector  with  integer  elements  , ,2 . 
Expanding  the  squared  refractive  index  n(x)  and  an  arbitrary  </>  in  Hilbert  space:  n(x^  - Yn  <p  (x)  fax)- Y&V  (x)’  where: 

reZ2  /eZ2 

\dxn(x)i=m(A)YJ\nl\1'  and  using  <PM)<PM)  = <P,aM)  for  Me  Z"  whereas,  Z 2 to  be  an  additive  group, 

A tzZ1 

fax)2  fax) = y[  Yn>-^s  ■ 

rezAseZ2  ) 


As  a result,  the  Helmholtz  (equation  (1))  with  ¥(x)  = for  some  given  k = ^ e B (with  T = (TiiT2) e z2 ) 

written  as  [2]: 

\7T2tT {Ar AY' t </),_,  = 1^,6,  + t2b2\\"  </>,_,  = 2 /7,|  ( 1 3) 

se/?2| 

Where:  te  Z2 ■ Equation  (13)  solved  for  {^)  in  complex  Hilbert  space  M(z2)  • 

This  system  written  as  below  [2]: 

...(14) 

4^r2 (r  + r)T (AT A)  '(t  + z)<p,  = ||(f,  + T,)b,  + (r2  + r,)bj[ </>,_,  = A Y ni+r-,A 

seR2 

Where:  m = A7r2diag  ([t  + rf  (AT A)_1[f  + r\)  , is  an  unbounded  diagonal  matrix  with  nonnegative  entries. 

T,  a multi-index  Toeplitz  matrix, <„(X)2  : xe  A | is  the  multi-index  spectrum.The  compact  subset  of  (0,+°o),  show  the 
condition ^(x)  as  a continuous  function  of  xs  A. 

Simulation  Result  and  Discussion 

This  simulation  results  design  ID  photonic  crystal  and  calculate  the  statistical  analysis  measurements  with  plots 
the  photonic  bands  for  2DPC.  Proceed  with  fields  for  TE  and  TM  modes.  The  results  below  are  got  after  following  these 
steps: 
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Found  the  total  number  of  plane  waves  used  in  Fourier  expansions. 

• Calculate  the  matrix  of  Fourier  coefficients  of  dielectric  function. 

• Calculate  the  points  along  the  perimeter  of  irreducible  Brillouin  zone. 

• Calculate  the  diagonal  matrices. 

• Solve  the  eigenvalue  problems  for  TEM  polarization. 

Figures  2 and  Figures  3 represent  the  digital  matrix  as  a wireframe  surface  for  TE  and  TM  polarization  photonic 
crystal  respectively.  The  TE  mode  statistical  analysis:  mean=0.4684,  median=0.4833,  mode=4.346E-005  and  the  standard 
deviation  (STD)=0.21 12.  Whereas,  mean=0.5218,  median=0.5425,  mode=4.953E-005  and  STD=0.2156  for  TM 
polarization  mode.  However  we  can  see  the  more  importanceof  polarization  effects  in  two  dimensions  systems  after 
compere  above  statistical. 

So  that,  in  order  to  understand  how  to  build  the  structure  of  the  optical  photonic  band;we  studied  a 2D  system.  Although 
the  sample  presented  in  this  paper  is  2D  photonic  crystal  is  easier  to  explain.  Calculations  were  limited  to  the  case  of 
Fourier  expansion  in  Hilbert  space  for  TE  and  TM  modes.  The  relationship  between  the  frequency  and  the  wave  vector 
(dispersion  relation)  in  vacuum  for  free  photons  is:  w=ck,  where  c is  the  light  speed  in  vacuum. 


(a) 


o.#. 

a.e. 


(b) 


Figures  2:  (a)  Display  a photonic  crystal  TE  polarization  digital  matrix  as  a wireframe  surface  for  three  variables 
matrices  defining  values  of  surface  crystal  components,  (b)  Display  a photonic  crystal  TM  polarization  digital 
matrix  as  a wireframe  surface  for  three  variables  matrices  defining  values  of  surface  crystal  components. 
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(b) 
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Figures  3:  (a)  Display  isolines  of  a photonic  crystal  TE  polarization.  A surface  represented  by  a digital  matrixplotted 
variables.  One  variable  ( coTE  )-plot  contour  lines  of  the  elements  of  a matrix,  interpreting  values  as  heights  with 
respect  to  the  x-y  plane.  Matrix  must  be  at  least  2-by-2  with  three  variables,  (b)  Display  isolines  of  a photonic 
crystal  TM  polarization.  A surface  represented  by  a digital  matrixplotted  variables.  One  variable  ( coTM  )-plot 
contour  lines  of  the  elements  of  a matrix,  interpreting  values  as  heights  with  respect  to  the  x-y  plane.  Matrix  must 

be  at  least  2-by-2  with  three  variables. 


The  modeled  structure  behind  Figures  4 and  Figures  5 are  a triangular  lattice  where  the  radius  is  r = 0.35.  This 
structure  presents  the  TE  and  TM  polarization  photonic  bands  for  2D  PC  consisting  of  cylinders  and  circular  cross-section 
and  infinite  height  arranged  in  a triangular  lattice.  The  corner  plots  represent  the  vectors  as  arrows  emanating  from  the 
origin  on  circular  grid  inside  the  photonic  crystal.  Two  variables  compass  graph  having  m arrows,  where  m is  the  number 
of  the  two  elements  variables.  The  tip  of  each  arrow  is  at  a point  relative  to  the  base.  This  syntax  is  equivalent  to  compass 
real  and  imaginary  elements. 


Figures  4:  The  TE  modes  photonic  crystal  band  structure.  The  right  corner  plots  represent  the  vectors  as  arrows 
emanating  from  the  origin  on  circular  grid  inside  the  TE  photonic  crystal. 


Figures  5:  The  TM  modes  photonic  crystal  band  structure.  The  right  corner  plots  represent  the  vectors  as  arrows 
emanating  from  the  origin  on  circular  grid  inside  the  TM  photonic  crystal. 


Impact  Factor  (JCC):  2.3529 


Index  Copernicus  Value  (ICY):  3.0 


Tem-Polarization  2D-Photonic  Crystal  Design  Using 
Fourieexpansion  in  Complex  Quantum  Hilbert  Space 


7 


CONCLUSIONS 

The  theoretical  modeling  and  simulations  are  presented.  At  the  beginning  the  background  theory,  which  is 
necessary  to  understand  photonic  crystals  presented.  Two  dimensional  photonic  crystals  are  investigated  by  using  both 
Fourier  expansion  and  Hilbert  space  method.  The  dispersion  relation  of  the  PCs  is  obtained  by  using  MATLAB. 

In  this  paper,  finding  a satisfactory  analytic almethod  for  the  spectrum  bandwith  the  mathematical  properties  of  two 
dimensions  optical  photonic  crystals  were  studied. The  structure  of  the  optical  band  gives  us  information  about  the 
electromagnetic  radiation  propagationproperties  inside  the  optical  photonic  crystal  and  when  the  plotof  energy  states  as  a 
function  of  propagation  direction  is  available,  the  structure  is  representation. 
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